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ABSTRACT 

A method for estimating redshifts of galaxy clusters based solely on resolved Sunyaev- 
Zel'dovich (SZ) images is proposed. Given a high resolution SZ cluster image (with FWHM 
of ~ 1'), the method indirectly measures its structure related parameters (amplitude, size, etc.) 
by fitting a model function to the higher order wavelet momenents of the cluster's SZ morphol- 
ogy. The applicability and accuracy of the wavelet method is assessed by applying it to maps 
of a set of clusters extracted from hydrodynamical simulations of cosmic structure formation. 
The parameters, derived by a fit to the spectrum of wavelet moments as a function of scale, are 
found to show a dependence on redshift z that is of the type x(z) = X\ exp(-z/x2) + JC3, where 
the monotony of this functional behaviour and the non-degeneracy of those parameters allow 
inversion and estimation of the redshift z. The average attainable accuracy in the z-estimation 
relative to 1 + z is ~ 4 - 5% out to z — 1.2, which is comparable to photometric redshifts. For 
single-frequency SZ interferometers, where the ambient fluctuating CMB is the main noise 
source, the accuracy of the method drops slightly to (A z/(l + z)) ~ 6 - 7%. 

Key words: galaxies: clusters: general, cosmology: distance scale, cosmology: cosmic mi- 
crowave background, methods: numerical 



1 INTRODUCTION 

Inverse Compton scattering of cosmic microwave background 
(CMB) photons off thermal electrons within the hot intra-cluster 
medium (ICM) of galaxy clusters produce fluctuations in the sur- 
face brightness of the CMB, an effect known as the thermal 
Sunyaev-Zel'dovich (SZ) effect (e.g. Sunyaev & Zel'dovich 1972, 
1980; Rephaeli 1995). Imaging clusters of galaxies through their 
SZ signature has, until recently, been a very challenging undertak- 
ing. To date, the development of detectors and new techniques have 
allowed high quality interferometric imaging of more than fifty 
clusters of galaxies (Carlstrom et al. 2002), despite incomplete cov- 
erage of the Fourier plane. In the foreseeable future, the availabil- 
ity of detectors in the microwave regime with angular resolutions 
surpassing 1' and sensitivities below jiK (e.g., the South Pole Tele- 
scope, described in detail in Carlstrom et al. 2002), will probe the 
hot plasma in galaxy clusters out to large redshifts providing SZ 
based wide field galaxy cluster catalogues and yielding a multitude 
of information about cluster formationand the cosmological model 
(Birkinshaw 1999). 

In particular, the abundance of clusters as a function of red- 
shift has been shown to be a very sensitive probe of the cosmologi- 
cal model (Eke et al. 1998; Henry 2000). The near independence of 
the line-of-sight SZ amplitude on cluster redshift makes the SZ ef- 
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feet the main tool for detecting galaxy clusters at high redshifts, 
0.5 i z ~ 2 (the upper limit depends on cosmology quite sensi- 
tively). This range of redshifts is especially important for probing 
the nature of the dark energy of the universe, since during this era 
it is expected to evolve rapidly until it eventually dominates over 
the other cosmological fluids. In order to obtain precise constraints 
on cosmological models it is essential to have accurate measure- 
ments of the redshift distribution of galaxy clusters (see Haiman 
etal. 2001). 

Normally, one determines the distance to the cluster by photo- 
metric or spectroscopic observations of the cluster member galax- 
ies. Unfortunately, this is a very challenging and time consuming 
task, in particular, when one considers the very large number of 
mostly high redshift clusters expected to be observed with sensi- 
tive future SZ instruments - The Planck satellite alone is expected 
to detect about 10 4 clusters (Bartelmann 2001). In order to replace 
photometric follow-ups we aim at inferring the distance to a clus- 
ter from SZ data alone for a future generation of experiments with 
increased angular resolution of about 1'. 

Theoretically, the cold dark matter (CDM) hierarchical clus- 
tering paradigm predicts a universal profile for dark matter ha- 
los that depends only on two parameters: core radius and density 
(Navarro et al. 1995). In addition, the same theory provides a very 
simple recipe for the mass accretion history of a certain halo as a 
function of its formation and observation redshift (Wechsler et al. 
2002; van den Bosch 2002; Zhao et al. 2003). Using these relations 
together with simple assumptions like hydrostatic equilibrium and 
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isothermality, one can expect that in the framework of the spheri- 
cal collapse model the observable SZ flux and apparent size should 
provide measures of the cluster's mass and distance. 

Indeed, using scaling relations, Diego et al. (2003) have 
demonstrated the viability of determining reliable morphological 
redshifts and examined different SZ observables with respect to 
their distance sensitivity. Among those observables, they showed 
that the cluster apparent size and central amplitude are promising 
distance indicators, once their degeneracy is broken. 

The main goal of this work is to derive redshifts of clusters 
based solely on their resolved SZ images by modeling the evolution 
of their structural parameters with redshift from the data set itself. 
This phenomenological approach does not depend on a priori as- 
sumptions about scaling relations that are valid only for spherically 
symmetric and relaxed systems. 

Specifically, the structural morphology of the cluster's pres- 
sure profile in an SZ observation is characterised by wavelet anal- 
ysis. 1 We are able to show that there is a simple relation between 
the distribution of moments over various scales in wavelet space 
and the cluster properties which can be described with simple 
phenomenological functions. Furthermore, the parameters of these 
functions are shown to follow a well defined and simple redshift de- 
pendence. Wavelet analysis has been chosen because it maintains 
the scale and positional information of cluster morphology, hence, 
it makes isolation and suppression of various unwanted contribu- 
tions to the observed signal possible while it reliabely upholds the 
underlying behavior. We note however, that Fourier space analysis 
could in principle yield very similar results. 

Hydrodynamically simulated clusters are used to demonstrate 
the method and to set limits on the redshift uncertainty expected 
in this approach. The simulated clusters used in the analysis are 
close to virialisation, e.g. merging systems are excluded. Under this 
restriction, both the relation between the observed quantity and the 
cluster physical parameters as well as the structural parameters are 
well defined. In addition, simulated clusters ignore radiative and 
feedback processes, the effect of which is discussed later in the 
paper. 

In the observational application, the evolution of the struc- 
tural parameters following from wavelet decomposition could be 
calibrated from a (relatively small) learning set of high quality SZ 
clusters with known (photometric/spectroscopic) redshifts. 

Our method relies crucially on the availability of resolved SZ 
cluster images. Therefore, throughout the paper we assume an in- 
strumental resolution of 1', where massive clusters should be re- 
solved even at the largest redshifts considered here. Indeed, future 
instruments such as the South Pole Telescope 2 (Carlstrom et al. 
2002) or the Atacama Cosmology Telescope 3 are designed to yield 
observations of up to 10 4 galaxy clusters with masses i 10 14 M Q 
(1 pK sensitivity) and x V resolution. 

This article is organised as follows: After basic definitions 
concerning the SZ effect in Sect. 2 and wavelets in Sect. 3, the 
simulations are outlined in Sect. 4. The capability of wavelets with 
respect to distance estimation is examined in Sect. 5. Possible sys- 
tematics are addressed in Sect. 6. A summary of the techniques in 
Sect. 7 and of the results in Sect. 8 concludes the article. 

1 There are also various ways of characterising the cluster's density profile 
in an SZ observation that are more or less susceptible to noise, for instance 
the fitting of a /3-prorile (Cavaliere & Fusco-Femiano 1978) to the electron 
density. 

2 http://astro.uchicago.edu/spt/ 

3 http : //www . hep . upenn . edu/ " angelica/act/act . html 



2 SUNYAEV-ZEL'DOVICH DEFINITIONS 

The SZ effect has been described in detail by many authors (for a 
comprehensive review see Birkinshaw 1999); here we briefly re- 
view its main aspects. The SZ effect arises because CMB pho- 
tons experience Compton-scattering off electrons of the dilute intra- 
cluster plasma. The CMB spectrum is modulated as photons are 
redistributed from the low-frequency part of the spectrum below 
218 GHz to higher frequencies. The change in thermodynamic 
CMB temperature due to the thermal SZ effect is 

Ar / e* + 1 \ 

— W)=yW)\ x — l -AU-2y(<l>) forx«l, (1) 

where x = hv/kT C MB is the dimensionless frequency. In the 
Rayleigh- Jeans limit (x <k 1), the change in temperature is asymp- 
totically equal to -2y($). The SZ amplitude at location 0, which 
is known as the Comptonisation parameter y(4>), is defined as the 
line-of-sight integral of the temperature-weighted thermal electron 
density: 

y(0)= f d/n e (<M)r e (<M). (2) 

where m c , c and k denote electron mass, speed of light and Boltz- 
mann's constant, respectively. r e (0, 1) and n e (0, 1) are electron tem- 
perature and electron number density at position (f> and distance /. 



3 WAVELETS 

3.1 Wavelet definitions 

During the last decade, wavelet analysis has become a popular tool 
in various data analysis and image processing applications. The 
main appeal of wavelet functional bases stems from their simulta- 
neous localisation of a signal in both the wavenumber and position 
domain, where they make orthogonal and complete projections on 
modes belonging to both spaces possible. In particular, the discrete 
wavelet families, by virtue of their constituting a complete basis, 
provide a unique and fast decomposition of the images into wavelet 
expansion coefficients. Statistics in terms of the q' h moments of 
the distribution of wavelet coefficients as a function of scale com- 
presses the signal contained in an image into a small number of 
parameters and yields information surpassing that derived in tradi- 
tional Fourier analysis. 

Following Daubechies & Bates (1993) and Muzy et al. (1993), 
the wavelet transform of a 2-dimensional image is defined as a 
convolution of the function y{x) to be analysed with the wavelet 
l/fr (\x - fi\): 

X Qi,o-) = J d 2 xy(x)- ^(\x-n\). (3) 

High values for x(t>, & ) are obtained in case of a match between 
the features of y(x) and the wavelet tfto-(x) at position fi and scale 
cr. From the wavelet expansion coefficients x(f*, cr) on scale cr at 
location fi one obtains the wavelet moments X q (cr) by integration 
over all displacements fi: 

X q (a) = J A 2 ii\x(n,o-)\". (4) 

The exponent q e N defines the order of the wavelet moment X q (cr). 
Values for q equal or larger than 2 allow noise suppression. The log- 
arithm In X(cr, q) of the wavelet moment as a function of logarith- 
mic scale In cr constitutes the wavelet spectrum. The A - 9 (cr)-statistic 
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is the main tool used in this study for characterising the morphol- 
ogy of SZ clusters. 

3.2 Application of wavelets to a cluster profile 

3.2.1 Analytic wavelet transform of a cluster y -profile 

In order to illustrate our idea of determining cluster sizes via 
wavelet decomposition, the wavelet transform of a King profile, 
which is known to describe the SZ morphology of clusters to first 
order, is performed. As an analysing wavelet, the Mexican-hat 
wavelet was chosen for simplicity. 

It is favourable to compute the convolution in the definition of 
X(ji, <r) in the Fourier domain. By virtue of eqn. (6), 

X(fi, cr) = J d 2 x y(x) \ff a (x - n) 

= {2nf J d 2 jfcy(*)Y -(-*;)exp(i*|i), 



(5) 
(6) 



the convolution reduces to a mere multiplication of the Fourier 
transforms Y(k) and ^(A:) of the image y(x) and the wavelet \l/ a {x), 
respectively. Restricting the order of the wavelet moment to q = 2 
and inserting the convolution theorem (6) into the definition (4) 
yields: 



X 2 (o-) = (2nf Jd 2 p d 2 kY(kyV a (-k)exp(ikn) 
= (2nf J d 2 k\Y(k)\ 2 \VAk)\ 2 , 



(7) 
(8) 



where the replacement ^(-A:)! 2 = ^(A)! 2 holds for real 
wavelets. 

The Mexican-hat wavelet is defined as the negative Laplacian 
of a Gaussian: 



^mhW = ^mh(*) 



1 



2mo- 2 6XP 2cr? 



(9) 



whereof the Fourier transform YmhW is derived by twofold partial 
integration: 

r d 2 x 

"^TTT '/'MhM SXp(-lkx) (10) 

J (2tt) 2 



(2. 



— f 



r dr(2o J2 -r 2 )exp\-—\J (kr) (11) 



(2n) 2 



■ exp 



k 2 cr 2 

2 



(12) 



where the azimuthal symmetry and the definition of the ze- 
roth order Bessel function of the first kind, 2nJo(kr) = 
j Q d(f> exp (ikr cos <f>) was used in the first step. Thus, the Fourier 
transform of the wavelet, 'i'a-(.k), is given by the Hankel transform 
of the Laplacian of a Gaussian. 

For the determination of Y(A), we assume that the projected 
thermal electron density can be described by a spherically symmet- 
ric King profile, i.e. a /3-model (Cavaliere & Fusco-Femiano 1978) 
with j3 = 1, core radius r c and central value of the Comptonisation 
parameter y : 



y(x) = y(r) = y Q 



1 + 



Then, the Fourier transform is easily computed: 
Y(k) = J ^y(*)exp(-iA*) 



(13) 



(14) 



y^i f 

2n J 



dr ■ 



ri + r 



Ukr) = ^ • K (kr c ), 



(15) 



where in eqn. (15) the definition of the zeroth order modified Bessel 
function of the second kind K (kr c ) was inserted. 

Substituting eqns. (12) and (15) into eqn. (8) and exploiting 
the azimuthal symmetry of the functions y(x) and ij/(x) yields an 
analytic integral for X 2 (o~): 



(16) 



X 2 (cr) = 2ny 2 r 4 c dk k 5 exp (-cr 2 A: 2 ) K 2 (kr c ). 

After evaluation of the integral in eqn. (16), the wavelet transform 
of the /3-profile reads as follows: 



X 2 (cr) 



2rl 



G xl (a 2 \- 2 'A 

^2,3 \ u I 000 J' 



(17) 



where a = r c jcr has been substituted. The function Q is Meijer's 
G-function, the exact definition of which is given by Gradshteyn 
& Ryzhik (1994). It is an interesting consistency to note that apart 
from the normalisation, the functional shape of eqn. (17) only de- 
pends on a, i.e. on the core radius r c expressed in units of the 
wavelet scale cr. 



3.2.2 Asymptotics of the analytical wavelet transform 

The asymptotic behavior of X 2 (cr) at the limit of cr <<c r c can be 
explored by substituting the expressions given in eqns. (13) and (9) 
into eqn. (6), and exchanging, by partial integration, the function 
on which the Laplacian operates. In the limit of interest the Gaus- 
sian can be replaced by a Dirac-(5 distribution. Substituting all of 
this into eqn. (4) yields that lim^o ^2(0") is proportional to y 2 Q and 
independent of cr, i.e. the normalisation of the wavelet spectrum 
measures the square of the central Comptonisation parameter y : 



Mo-) 



32jt y 2 , 
15 ' rl 



for 



a <K r c . 



(18) 



In the opposite limit, i.e. r c <<c cr, one can use the fact that the King- 
profile is highly peaked at the center and that it is convolved with a 
Mexican-hat wavelet guaranteeing the convergence of the integral 
in eqn. (16) at 00. In the limit of r c — > this integral is dominated 
by the value at k = 0. Therefore, one can approximate the King- 
profile with a Dirac-(5 distribution and show that asymptotically the 
lim^o X 2 (c) is proportional to cr' 6 : 



X 2 (cr) cc 



for 



cr » r c . 



(19) 



The sensitivity of the wavelet spectrum X2(cr) on cluster size 
r c is illustrated by Fig. 6. The wavelet spectrum is constant for cr <<c 
r c , has an r c -dependent break and drops off asymptotically cc cr' 6 
for cr » r c . Naturally, the scale cr at which the transition from one 
asymptotic regime to the other occurs, is determined by the value 
of r c , i.e. the cluster size. 

Motivated by this example, the wavelet moments X q {cr) ob- 
tained from real data (Sect. 5.3) will be fitted with a power law 
with an exponential cutoff, where the cutoff indicates the cluster 
size and the amplitude is proportional to some power of the central 
Comptonisation parameter y . 



3.2.3 Finite instrumental resolution 

The influence of finite instrumental resolution can easily be incor- 
porated by an additional factor \B(k)\ 2 in eqn. (8): 



4 B. M. Schdfer, C. Pfrommer and S. Zaroubi 




wavelet scale cr 



Figure 1. Sensitivity of the wavelet spectrum on the cluster size: The second 
order wavelet moments X2(<r) are shown as a function of cr for various core 
sizes r c = 0.5 (solid), r c = 1 (dashed), r c — 2 (dotted) and r c = 4 (dash- 
dotted). The curves have been normalised to their asymptotic values for 

(T -> 0. 



X 2 {cr) = {2nf J d 2 k\Y(kf \^Ak)\ 2 \B(kf , 



(20) 



where B(k) is the Fourier transform of the (azimuthally symmetric) 
beam profile b(x), which is for simplicity assumed to be of Gaus- 
sian shape with FWHM = V8 hi(2) • a b : 

B(k) = J ^&(x)exp(-ifa;)with (21) 



b(x) 



x 



(22) 



* \ b ' 

This effectively replaces cr in eqn. (17) by the harmonic mean 

yjcr 2 + cr\, which limits the range of accessible wavelet scales to 
cr > cr b . 



3.3 Analogy to power spectra in Fourier analysis 

By interpreting the wavelet spectrum in eqn. (8) as the variance of 
the fluctuations on the scale cr, one may draw an analogy to Fourier 
decomposition: 



(23) 



var \y(x)] = X 2 (cr) = (2nf J d 2 k P(k) |>P 

where P(k) = (27r) 2 (\Y(k)\ 2 ) is the Fourier power spectrum. The 
wavelet tf/(x) now adopts the role of a filter function on scale cr. 
This filter function reads in real space in the case of the Mexican 
hat wavelet: 



^mhO*) = 



2^-cr 6 



■ ■ exp 



2cr 2 ) 



(24) 



Therefore, our method is equivalent to considering power spectral 
analysis of filtered fields and higher order Fourier space moments. 



4 SIMULATIONS 

The accuracy in the determination of redshift z was assessed by 
examining the performance on numerical simulations: First, simu- 
lations of cosmological structure formation including gas physics 



have been carried out in order to model the evolution of clusters 
(Sect. 4.1). Subsequently, maps of the Compton-y parameter have 
been produced by using an interpolation kernel with an adaptive 
smoothing length for projecting the Compton-y parameter along 
the line-of-sight (Sect. 4.2). By applying selection criteria favour- 
ing virialised systems a cluster sample was compiled (Sect. 4.3). 
Finally, aiming at realistic single frequency SZ observations, we 
simulated the ambient CMB fluctuations that act as the primary 
source of noise (Sect. 4.4) and combined the resulting realisations 
of the CMB with the cluster maps (Sect. 4.5). 

The assumed cosmological model is the standard ACDM cos- 
mology, which has recently been supported by the WMAP satellite 
(Bennett et al. 2003; Spergel et al. 2003). Parameter values have 
been chosen as fi M = 0.3, fi A = 0.7, H = 100/ikms" 1 Mpc~' 
with h = 0.7, n B = 0.04, n s = 1 and cr 8 = 0.9. 

4.1 SPH cluster simulations 

A simulation of cosmological structure formation kindly provided 
by V. Springel and L. Hernquist (Springel & Hernquist 2002; White 
et al. 2002) constitutes the basis of our analysis. In a cubical box 
of comoving side length 100 Mpc/h with periodic boundary condi- 
tions a smoothed particle hydrodynamic (SPH) simulation compris- 
ing 216 3 dark matter particles as well as 216 3 gas particles was run 
and snapshots were saved at 23 redshifts ranging from z = 0.102 
out to z = 1.114. The comoving spacing along the line-of-sight of 
two subsequent boxes is 100 Mpc/h. Purely adiabatic gas physics 
and shock heating were included, but radiative cooling and star for- 
mation were ignored, which however does not result in significant 
differences in SZ morphology, as shown by White et al. (2002) but 
impacts on the scaling relations as demonstrated by da Silva et al. 
(2001). 

Overdensities are identified using a friends-of-friends algo- 
rithm with the linking length b = 0.164, which yields all member 
particles of a cluster in conjunction with a spherical overdensity 
code, from which virial quantities are estimated. We computed the 
mass M vir inside a sphere of radius r vir , interior to which the aver- 
age density was 200 times the critical density p crit = 3H 2 /(8nG). 
The angle subtended by twice the virial radius is denoted as £> vir . 
We imposed a lower mass threshold of M vir > 5 • 10 13 Mq/H. 



4.2 SZ map preparation 

Square maps of the Compton-y parameter of the selected clusters 
were generated by SPH projection of all member gas particles onto 
a cubical grid with 128 2 mesh points. The (comoving) side length 
.v of the map was adapted to the cluster size, such that the comov- 
ing resolution g = .v/128 of the grid is specific to a given map. 
Examples of Sunyaev-Zel'dovich maps are given in Fig. 2. 

If the particle p at position r p = (x p ,y p ,z p ^ has a smoothing 
length hp, an SPH electron number density estimate n p , and an SPH 
electron temperature T p , the Compton-y parameter at the pixel at 
position x is given by: 



= — 2~r / i 



x+g/2 y+g/2 



Jdx p Jdy p jdz p w {j~^ ' "p T p 



x-g/2 y-g/2 



with r = ^Jix,, - xf + (y p - y) 2 + zj. 



(25) 



(26) 



Here, we assumed complete ionisation and primordial element 
composition of the ICM for the determination of electron number 
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Figure 2. Picture book of Sunyaev-Zel'dovich clusters: The upper panel shows clusters at high redshifts of z — 0.8, in comparison to clusters at low redshifts of 
Z = 0.174 in the lower panel. The columns contrast different morphologies in an exemplary fashion: relaxed systems (left column), elongated clusters (centre 
column) and clusters in the phase of minor merging or mass accretion (right column). The grey scale denotes the amplitude of y(x)' = log [ 1 + 10 5 ■ y(jc)| and 
the contours have a logarithmically equidistant spacing of 0.1 dex, i.e. the lowest contour denotes a common value of y = 2.5 ■ 10~ 6 . 



density and temperature. In this way we produced projections along 
each of the three coordinate axes. The function w is the spherically 
symmetric cubic spline kernel suggested by Monaghan & Lattanzio 
(1985), which was also used in the SPH simulation: 

J" 1-6m 2 + 6m 3 ,0<m< 1/2 
w(u) = --\ 2(1 -uf ,l/2<u<l with u = f/hp. (27) 
n \ ,u> 1 

The fact that the kernel w is defined on a compact support 
u G [0 ... 1] greatly reduces the computational effort. 

4.3 Cluster selection 

Clearly, the wavelet redshift estimation relies on the clusters not 
being in the state of violent merging. Apart from the minimal mass 
of M min = 5 ■ 10 13 Mq//!, that translates into a minimally required 
line-of-sight Comptonisation amplitude y mm , clusters have been se- 
lected in order to meet the following prerequisites: 

• By visual inspection it was made sure that the clusters show 
but a single peak in the Compton-y map in order to exclude systems 
in the late phase of a merger. 

• The SZ morphology is required not to be too elongated. By fit- 
ting a 2-dimensional /3-model yp{x) to the SZ profile yamaC*), values 
for the semi-axes r x and r y are derived. 90% of the clusters within 
the selected sample have axis ratios q = r y /r x smaller than 0.8 and 

ellipticities e = ^jr 2 x - r^/r x below 0.6. 

• Residual deviations from the canonical /3-profile ought to be 



small. The rm.s-deviation v of the cluster from the best-fitting /3- 
profile, 

was smaller than 25% for 90% of our clusters. 

Applying these selection criteria, 10 clusters were selected 
from each of the 23 redshift bins, yielding with the three orthog- 
onal projections of each cluster a total number of 690 maps with 
which the accuracy of the wavelet method in estimating redshifts 
was assessed. The distributions of the ellipticities e and the inte- 
grated residuals v are shown in Fig. 3. The same distributions were 
derived for the smoothed cluster maps, were the effects of finite in- 
strumental resolution have been incorporated. As Fig. 3 suggests, 
the beam does not have a major impact on the morphological prop- 
erties of the cluster sample, which is due to its narrowness of only 
1' (FWHM). 

4.4 CMB map generation 

CMB anisotropies are assumed to be a particular realisation of a 
Gaussian random field. Aiming at simulating a realisation of the 
CMB on a square, flat map, we take temperature fluctuations (<f>) 
relative to the average CMB temperature of (T) = 2.726K to be the 
independent random field, 

T(<b)-(T) 

9(4>)= {T) ■ (29) 
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Figure 3. Selection criteria: distribution of residual deviations v from the 
best-fitting yS-profile for unsmoothed (solid, circles) and smoothed (dashed, 
diamonds) maps. The second set of lines shows the distribution of the ellip- 
ticies e with (dash-dotted, squares) and without smoothing (dotted, crosses). 

The flat, two-dimensional power spectrum P g (£) is defined via: 

(&(t)@'(e')) ee (2nr 2 s D ({ - e')p e (\e\), (30) 

where &(f) denotes the Fourier transform of 6($). The simulation 
of the CMB temperature fluctuations on a flat square map now con- 
sists of the following two steps: 

• The angular power spectrum Ce is computed for the flat 
ACDM-Universe using the CMBfast code by Seljak & Zaldarriaga 
(1996). In addition to the cosmological parameters being already 
described in Sect. 4, we use adiabatic initial conditions and set the 
primordial He-mass fraction to X He = 0.24 and the Thomson op- 
tical depth to r = 0.17 (Spergel et al. 2003). The angular power 
spectrum of the CMB is normalised to COBE data. Since the SZ 
effect distorts the CMB only on small angular scales, the flat sky 
approximation I » 1 is fulfilled and it is appropriate to replace 
the spherical harmonics with plane waves. Hu (2000) has shown 
that the 2-dimensional flat power spectrum P e (i) is approximately 
equal to its angular analogue: Ct - Pgit). 

• Then, Gaussian random variables are generated on a com- 
plex two-dimensional grid in Fourier space with variance cr 2 (() = 
Pe(0 according to the absolute value of their wave vectors t. 
Inverse Fourier transform yields a realisation of the temperature 
anisotropies 6(<f>). 

4.5 Simulated single-frequency SZ observations 

For SZ clusters observed with a single-frequency interferometer 
(e.g., the CBI experiment, Halverson et al. 2002) 1 , it is important to 
examine the applicability of the X 9 (<x)-statistic. For the purpose of 
this paper, it suffices to consider observations at small frequencies 
v. Thus, the Compton-y maps are combined with realisations of the 
CMB fluctuations by taking advantage of eqn. (1) in the Rayleigh- 
Jeans limit, 

7/(0) = [l-2y(0)][l+ 0(0)] <r>. (31) 
1 http : //www . astro . caltech . edu/~ t jp/CBI/ 
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Figure 4. Simulated temperature map of the CMB combined with a fore- 
ground SZ cluster atz = 0.102 with virial quantities M vir = 2.2- lO I4 M0//i, 
r v j r = 1.47 Mpc//i and fc7\,j r = 1.52 keV. At the cluster centre, the SZ tem- 
perature decrement amounts to -1.8 mK and the CMB temperature fluctu- 
ation with the highest amplitude is equal to 0.23 mK. A total of 30 linearly 
spaced isothermals are drawn. In this case, the comoving scale 1 Mpc//i 
corresponds to 11.5'. 



Fig. 4 shows the Compton-y map of a nearby cluster of 
2.2 • \0 l4 M Q /h at redshift z = 0.102 combined with a patch of 
the CMB constructed by the algorithm described above. In this 
map, the average CMB temperature (77) was subtracted. In order 
tomimic observations, the resulting combined maps are smoothed 
with a Gaussian beam with FWHM of V8 ln(2) ■ cr b = V. 

In the case of multi-frequency SZ observations the SZ signa- 
ture can be easily distinguished from the CMB signal. Therefore, 
for these cases the CMB background is ignored and not included in 
the simulated cluster SZ images. Nevertheless, finite instrumental 
resolution was taken care of and the SZ maps were convolved with 
a Gaussian kernel of V8 ln(2) • cr b = V (FWHM). 



5 ANALYSIS 

In this section, the analysis is explained step by step: After in- 
troducing the wavelet families (Sect. 5.1), the wavelet spectrum 
and the parameters deduced from it are described (Sect. 5.2 and 
Sect. 5.3). The correlations of the wavelet spectral parameters with 
physical quantities are discussed (Sect. 5.4). The measurement 
principle and the breaking of degeneracy is illustrated in Sect. 5.5. 
Successively, the intercorrelation of the wavelet parameters and the 
shape of the parameter space is explored by principal component 
analysis (Sect. 5.6). Then, gauge functions for modelling the red- 
shift dependence ofthe parameters are proposed (Sect. 5.7). An im- 
portant issue for single-frequency interferometers is the influence 
of primary CMB fluctuations on the wavelet spectrum and their 
suppression (Sect. 5.8). Finally, the redshift of the clusters are esti- 
mated by maximum likelihood techniques (Sect. 5.9). 
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Figure 5. The wavelet basis functions ip(x) chosen for the analysis: symlet 
sym2 (solid), Daubechies' wavelet db4 (dashed), coiflet coifl (dotted) and 
the biorthogonal wavelet biorl.3 (dash-dotted). 



5.1 Wavelet basis functions 

In the analysis a wide range of wavelets with different functional 
shapes was employed, although the symlet wavelet basis introduced 
by Daubechies & Bates (1993) yielded particularly good results. 
Due to their symmetry and peakiness, symlets are seemingly es- 
pecially suited for analysing SZ morphologies, because spiky fea- 
tures are useful for edge detection, and since apparent cluster size is 
the primary observable, the usage of a peaked wavelet is justified. 
Other wavelet families that found application in our analysis were 
Daubechies' wavelets, coiflets and biorthogonal wavelets. Fig. 5 
compares the functional shape of the different wavelet families. 

The analysis proceeds by measuring the wavelet moments on 
smoothed comoving maps of the Compton-y parameter following 
the definition in Sect. 3. The scale cr of the resulting wavelet spec- 
trum is then converted to angular units. Because our SZ maps are 
computed on a grid of 128 2 mesh points with adaptively chosen side 
length for each cluster, our dynamical range of the wavelet spectra 
always comprises approximately two decades. However, this is no 
fundamental limitation of this approach because the maps are fea- 
tureless below the smoothing scale of 1' (FWHM). 



5.2 Measurement of wavelet quantities 

In order to derive the actual flux decrement or, equivalently, the 
decrement in antenna temperature from the line-of-sight Compton- 
y amplitude, its value at each pixel needs to be multiplied with 
the solid angle it subtends. For the conversion, a standard ACDM- 
cosmology was assumed, the parameters of which have already 
been given in Sect. 4. Thus, the pixel amplitudes were modified 
according to: 



y(x) — > y((/>) = y(x) ■ 4 arctan 



2w(z) 



(32) 



where w(z) is the comoving distance in the model cosmology and 
g denotes the comoving size of a single pixel. It should be empha- 
sised that the wavelet coefficients x(ft, cr) are evaluated on a comov- 
ing grid, which has been adapted to the cluster size before convert- 
ing the wavelet scale cr to angular units. This, however, should not 
pose a problem for real observations, provided the sampling scale 
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Figure 6. The spectrum of wavelet moments, together with the fitting for- 
mula (33) for increasing wavelet moment order q: q = 2 (squares), q = 3 
(circles, solid), q = 4 (stars, dashed), q = 5 (diamonds, dotted) and q = 6 
(crosses, dash-dotted) for a single cluster. The wavelet moments X q (cr) fol- 
lowed from wavelet expansion with the ,sywi2-wavelet. 



is of the same order of magnitude as the angular scale of our finest 
pixels. 

In order to obtain dimensionless quantities, the unit of the 
wavelet (fro-(x) has been set to inverse steradians, such that the 
wavelet expansion coefficients cr) and the wavelet moments 
X q {cr) are dimensionless, irrespective of q. For numerical conve- 
nience, the pixel amplitudes in the combined SZ maps have been 
multiplied with 10 12 . 

The summation in the definition of the wavelet moment X q (cr) 
in eqn. (4) discards the information about the position fi at which 
the wavelet expansion coefficient x(M>°~) ' s evaluated. Conse- 
quently, the position of a cluster inside the observing frame does 
not influence the wavelet decomposition. 



5.3 Wavelet spectrum of SZ cluster maps 

Due to the lack of any analytical generalisation of eqn. (17) for 
q + 2, deviations of the Compton-y map from a King profile and 
wavelets other than the simple Mexican hat, we decided to explore 
phenomenological functions for describing the wavelet spectrum. 
The simplicity of the shape of the wavelet spectrum shown in Fig. 6 
implies that the model function, 



lnX q (a) =s a + s In (tr/tr ) - cr/c, 



(33) 



is able to extract all apparently contained information, i.e. the spec- 
trum is described by means of three quantities: the amplitude a, the 
slope s and a break at c. The parameter cr has been included in 
eqn. (33) in order to obtain a formula that is dimensionally correct, 
although it does not yield any new information and this specific 
degree of freedom is already described by the variable a. 

The usage of eqn. (33) implicitly discards information about 
asphericity and effectively determines an average of the cluster's 
extension along its major axes. The problem would be significantly 
complicated by including asymmetry and considering vectorial na- 
ture of cr (see Zaroubi et al. 1998, 2001). 

Because the cutoff parameter c is of great importance to our 
analysis, it needs to be derived reliably. Thus, the order of wavelet 
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Figure 7. Wavelet measured cluster size c versus angular extension VIC for 
increasing wavelet moment order q: q = 3 (circles) and q = 6 (crosses) 
without including CMB fluctuations. The c-values have been determined- 
with the ,sy»i2-wavelet. 
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Figure 8. Wavelet amplitude a as a function of integrated Comptonisation 
parameter T for different weighting exponents q: q = 3 (circles) and q = 
6 (crosses), again without taking CMB fluctuations into account. As the 
analysing wavelet, the .vym2-wavelet was chosen. 



moments q was restricted to q > 3, because larger ^-values facili- 
tate the determination of c. From Fig. 6 it is obvious that an increase 
in q suppresses the value of X q (cr) at small scales cr such that the 
curve develops a maximum in the vicinity of c. Additionally, by 
the choice of large values for q, the wavelet expansion coefficients 
Xifi, cr) dominated by CMB noise are suppressed relative to those 
obtained in the central part of the cluster and consequently higher 
order wavelet moments X q (cr) provide a cleaner measurement. The 
range of sensible ^-values is restricted by the fact that for increas- 
ing q the moment X q (cr) is successively dominated by the largest 
wavelet expansion coefficient x(n, cr) and does no longer contain 
information of the structure to be analysed. In order to stabilise the 
fitting procedure we interpolate in between the wavelet moments 
X q (cr). This is justified because we expect a smooth variation of the 
wavelet spectrum according to Sect. 3.2.1. 



5.4 Correlations with physical quantities 

The parameters derived from the fit to the spectrum of wavelet coef- 
ficients have a physical interpretation. As illustrated in Sect. 3.2.1, 
the wavelet spectrum breaks at the cluster scale. 

Therefore, one expects a correlation between the angular size 
of the cluster B vir and the cutoff c, as shown by Fig. 7. Increasing 
weighting exponents q shift the regression line to smaller values of 
c, which can be understood by the fact that larger values of q sup- 
press small wavelet expansion coefficients arising at the outskirts of 
the cluster, which in turn leads to a break in the wavelet spectrum 
at smaller scales. 

Similarly, the amplitude a determined by the fit is propor- 
tional to the integrated Compton-y flux, 



T= d 2 <M0): 



kT viT a T l+ fn j M vi 
m e c 2 ^a(z) 2 2 m„ 



(34) 



as illustrated by Fig. 8. Here, f b denotes the baryon fraction, f H 
the hydrogen fraction, which determines the elemental composi- 
tion and has been set to the primordial value of 0.76, and m p is the 
proton mass. d A (z) is the angular diameter distance in our cosmol- 
ogy- 



The normalisation a of the wavelet moments X q (cr) shows a 
steeper dependence on the integrated Comptonisation parameter T 
for larger choices of q, which is explained by the following argu- 
ment: The amplitude a(q) reflects the normalisation of the wavelet 
moments X q (cr). The integral in eqn. (4) is dominated by the largest 
wavelet expansion coefficient cr), taken to the q lh power. On the 
other hand, the wavelet expansion coefficients xif 1 , cr) are propor- 
tional to integrated Comptonisation Y, resulting in observed rela- 
tion In [^(cr)] oc a cc q ■ ln(T). 

The slope 5 is a measure of instrumental smoothing: Placing 
the same cluster at different redshifts would result in a blurred im- 
age of the more distant one. Keeping in mind that there is a close 
analogy between wavelet- and Fourier-transforms (as explained in 
Sect. 3.3), the wavelet moment X q (cr) as a function of cr can be 
interpreted as the variance of the wavelet-filtered field. The instru- 
mental beam introduces an additional filtering to the Compton-y 
map (compare Sect. 3.2.3) and would cause the Fourier spectrum 
to drop at smaller values of the wave vector k, because the instru- 
mental beam constitutes effectively a low-pass filter that is erasing 
structures smaller than its extension. Comparing clusters at differ- 
ent redshifts, it is clear that the drop in power happens at smaller 
scales in the case of the more distant cluster. Then the slope s, de- 
fined as d In X q (cr)/ diner for cr «: r c , is larger in the case of a 
resolved cluster compared to an unresolved cluster. This measure 
of the influence of finite instrumental smoothing varies only by a 
factor of two in the redshift and mass range considered here, but 
nevertheless serves as an indicator of cluster distance. 



5.5 Measurement principle 

Now, it is necessary to illustrate how a measurement of the total 
Comptonisation Y and of the angular size 9 vlI suffices to derive a 
distance estimate. For that purpose, clusters are placed at unit dis- 
tance and the distance dependences of the wavelet amplitude a and 
the the cutoff c are removed by the following formulae, since a is a 
logarithmic measure of flux inside an solid angle element Y and c 
is a logarithmic measure of of angular extension vir : 



a = a(z) + 2 -In (d a (z)) 



(35) 
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Figure 9. Distance corrected wavelet amplitude ao = a{z) + 2 ■ In (d a (z)) as a 
function of distance corrected wavelet cutoff parameter cq = c(z)+ln (d a (zj). 
The values have been determined in fits to the wavelet spectrum X^cr), that 
has been derived with the sym2-wavelet as the analysing wavelet. 

co = c(z) + In (d a (z)). (36) 

By applying simple scaling arguments, one expects the ratio 
ao/t'o to be equal to 5: From the wavelet amplitude a one obtains 
a oc In (T • d A (z) 2 ^ oc In (M vir • r vir ). Furthermore, from the spheri- 
cal collapse model follows, that oc M 2 ' 3 (Navarro et al. 1995), 
which yields, together with M vir cc r 3 vir , the relation a oc ln(r^ ir J. 
Substituting Co oc In (r vir ) gives the final result a /c = 5. 

Fig. 9 nicely illustrates how the degeneracy is broken and how 
a simple measurement of flux and angular extension suffices in or- 
der to derive a distance estimate: A crude fit to the distance cor- 
rected wavelet amplitude a as a function of distance corrected 
wavelet cutoff parameter c yields a slope of approximately 5.8, 
which corresponds well to the slope of ~ 5 expected from the 
theoretical consideration outlined above. If, hypothetically, the ra- 
tio ao/ c 'o was equal to 2, the measurements of flux and angular 
size would be completely degenerate and would not yield any dis- 
tance information. This case corresponds to disks of equal surface 
brightness, where measurements of flux and angular size are com- 
pletely degenerate and do not yield any distance information at all. 
It should be noted, that by adopting the usual scaling relations, one 
introduces a systematic error in slope that can amount to ^ 20% 
error. 



5.6 Principal component analysis 

In order to rate the extent to which the parameters derived in the 
fit to the spectrum of wavelet moments are independent, a principal 
component analysis (PC A, see e.g. Deeming 1964) was performed. 
The PCA is determining a transformation to a new orthogonal co- 
ordinate system in parameter space spanned by a, c and s, such that 
the variance along the first axis is maximised. 

The first eigenvector of the matrix that describes the change 
of basis by the PCA reads x PCA = (0.65, 0.70, 0.32) which has been 
derived for the spectral parameters for q = 3 and with the sym2- 
wavelet as the analysing wavelet. The values similar in magnitude 
suggest that the variation in the data set is contained in all three 
parameters a, c and s. 



q = 3 q = 4 q = 5 q = 6 

1 st principal component 95.6% 94.2% 92.8% 91.5% 
2 nd principal component 2.7% 4.2% 5.5% 6.7% 



Table 1. Result of the PCA. The variance explained by the first and sec- 
ond principal component as a function of wavelet order q. Here, no CMB 
fluctuations were included. 

As can be read off from table 1, the parameter space is tightly 
constrained and all three parameters are interrelated, such that the 
data points form a narrow ray in parameter space. This result holds 
irrespective of the choice of q, although the scatter increases with 
higher choices for q. The values in table 1 have been determined 
without considering CMB fluctuations. Given the physical inter- 
pretation of the wavelet amplitude a and the cutoff c, it is obvious 
that the tight correlation can be traced back to the self-similarity of 
clusters and the cluster scaling relations linking r vir , M vir and r vir 
that follow from the spherical collapse model. The scaling relations 
for SZ quantities derived by da Silva et al. (2003) support this view. 
This shows together with Sect. 5.4 and Sect. 5.5 hat both the cutoff 
c and the amplitude a are functions of M vir . By only considering 
systems close to virialisation we thus expect a fundamental narrow 
ray in parameter space which can be described by a single principal 
component. 

5.7 Redshift dependence of the wavelet parameters 

The parameters a and c are expected to decrease with increasing 
redshift z, the reason for which is quite apparent: The angular di- 
ameter 6^ and the integrated Comptonisation T decreases because 
of the increasing angular diameter distance d A (z) that enters vir lin- 
early and T quadratically. Furthermore, clusters accrete matter dur- 
ing their formation history and thus are on average more massive at 
later times, i.e. at smaller redshifts z (see, e.g., Wechsler et al. 2002; 
van den Bosch 2002; Zhao et al. 2003). From the physical point of 
view, the dependence of a and c on redshift z is far from trivial, and 
therefore, their functional behaviour is described by an empirical 
approach. Among others, the exponential function provides a good 
fit to the data, as illustrated by Figs. 10 and 11: 

x(z) = x\ exp I - — | + x 3 , where x e {a, c, s] . (37) 

The optimised parameters xu ; e ( 1, 2, 3), for x e {a, c, s] in the 
gauge functions eqn. (37) are given in table 2 for the case q = 3. 
It should be emphasised that the parameters stated are only valid 
for image analysis with the .s>m2-wavelet, where the maps have 
been smoothed with a Gaussian kernel with 1' (FWHM) and the 
considered cluster sample, which is defined by the selection criteria 
laid down in Sect. 4.3 and the minimal mass of 5 • 10 13 Mq/A. 

5.8 Influence of CMB noise 

Clusters at high redshift z are characterised by their small angular 
scale on which the underlying CMB is represented by a smooth 
gradient due to Silk damping (Silk 1968). In this case the wavelet 
analysis produces the same results irrespective of the CMB noise 
owing to the distinct morphological feature of the cluster on top of 
the smooth CMB gradient. Once clusters at lower redshifts reach 
angular sizes comparable to characteristic scales of CMB fluctua- 
tions, the wavelet analysis has to be made more sophisticated. This 
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Figure 10. Dependence of the wavelet cutoff parameter c on redshift z with- 
out considering CMB fluctuations for q = 3 (circles, solid) and q = 6 
(crosses, dash-dotted). The analysing wavelet was the .vym2-wavelet. 



16 




-2 - 



-4 ' 1 ' 1 1 1 

0.2 0.4 0.6 0.8 1 1.2 

redshift z 

Figure 11. Dependence of the wavelet amplitude parameter a on redshift 
Z without including CMB fluctuations for q = 3 (circles, solid) and q = 6 
(crosses, dash-dotted). The quantities have been determined with the sym2- 
wavelet. 



parameter 


variable 


i= 1 


i = 2 


i = 3 


amplitude 


a 


10.5837 


0.6475 


-1.9570 


cutoff [arcmin] 


c 


0.5124 


0.5165 


0.3809 


slope 


s 


1.3423 


0.4144 


1.3803 



Table 2. Fitting values for the gauge functions defined in eqn. (37) for the 
cluster sample at hand and the .vym2-wavelet basis. The order of the wavelet 
moment X q (<r) has been set to q = 3. The values have been derived without 
taking CMB fluctuations into account. 
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Figure 12. Changes to the wavelet spectrum of a single cluster caused by the 
fluctuating CMB: unperturbed wavelet spectrum of the SZ cluster (circles, 
solid), of the pure CMB (squares, dash-dotted) and of the combined map 
(diamonds, dashed). Data points were derived from simulated data and the 
joining line in the case of the unperturbed wavelet spectrum the result of the 
fitting functions described by eqn. (33). The order of the wavelet moment 
is q = 6 (thick) and q = 4 (thin). Again, the analysing wavelet is the sym2- 
wavelet. 



complication in the wavelet analysis arises because wavelets are 
primarily suited for determining morphological features rather than 
solely singling out high amplitude characteristics. Because the an- 
gular scale of the clusters ranges between 10' and 1', which cor- 
responds to multipole orders of £ ^ 10 3 . . . 10 4 , it suffices to con- 
sider the Silk damping tail of the angular power spectrum of the 
CMB. In the wavelet spectrum X q (cr) this translates into an addi- 
tional approximate power-law component X q MB (cr), as can be seen 
from Fig. 12: 

In X™ B (o-) - flcMB + scmb In a. (38) 

This is due to the discrete sampling of the wavelet moments 
as well as the inherent statistics of the wavelet spectra of order q 
which can be interpreted as suitably weighted g-point correlation 
functions in Fourier analysis (compare Sect. 3.3). 

Fig. 13 shows the probability distribution function 
p(scMB)d.?cMB of the slopes .s C mb following from linear fits 
to the wavelet moments X q (a) for the range of q's considered here. 
Again, the .sym2-wavelet was chosen as the analysing wavelet. 
The slopes i C MB are not well confined, keeping the vast range 
of angular scales in mind, which in turn will make it difficult to 
subtract the CMB -contribution to the wavelet spectrum of the 
combined map. 

In order to disentangle the contributions from the CMB noise 
from those of the cluster, one may pursue different approaches: 
Among others, CMB fluctuations underneath the cluster can be re- 
constructed with spline polynomials and successively subtracted. 
Here, we have masked the cluster and fitted 5 th -order polynomi- 
als to the remaining data points. Because the y-maps and the re- 
alisations of the CMB are to leading order combined linearly and 
because the CMB is a smoothly varying field, it is possible to re- 
construct the CMB fluctuations from the environment of the cluster 
and interpolate to the cluster centre. The reconstructed CMB field 
can be subtracted from the inital image and by applying wavelet 
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Figure 13. Distributions of the power-law slopes .S'cmb of CMB wavelet 
spectra as a function of wavelet moment order: q = 3 (circles, solid), q = 4 
(stars, dashed), q = 5 (diamonds, dotted), q = 6 (crosses, dash-dotted). For 
the analysing wavelet, the symlet sym.2 was chosen. 

decomposition to the cleaned field one obtains a wavelet spectrum, 
from which the parameters a, c and s can be reliably derived. 

5.9 Redshift estimation 

In order to assess the accuracy of the redshift measurement, a max- 
imum likelihood estimation is performed. The likelihood function 
is defined as: 

L{Z) = (2n)^a a a c a, 6XP " P N If £H ' (39> 

and was evaluated for each bin separately, i.e. the index ( enumer- 
ates clusters within the redshift bin under consideration. N = 30 
denotes the number of clusters within a single redshift bin. From 
the position of the maximum in £(z) the most probable redshift es- 
timate z was derived and the accuracy of the estimate followed from 
the percentiles corresponding to 1-tr confidence intervals. Fig. 14 
shows the estimated redshift versus the real redshift for the cluster 
sample derived by using all of the three parameters a, c and s. In 
comparision, the error bars have become larger by a factor of ^ 1.5 
when including the fluctuating CMB, as illustrated by Fig. 15. The 
measurement is unbiased and the error relative to 1 +z rises slightly 
with increasing redshift z. 

The results for different analysing wavelets as a function of 
wavelet moment order q are summarised in tables 3 and 4. Clearly, 
the method starts to fail at redshifts exceeding i 1, when the angu- 
lar diameter distance dpSf) develops a plateau and does not cause 
clusters to appear smaller. The average attainable accuracy is stated 
relative to 1 + z in order to facilitate comparison to photometric 
redshifts. The accuracy slightly degrades with increasing q, which 
is due to suppression of small wavelet expansion coefficients es- 
pecially at small scales and the resulting inaccuracy of the fitting 
formula eqn. (33) used to extract the spectral parameters a, c and s 
from the wavelet spectrum. 

Inclusion of the CMB in order to test the applicability of de- 
termining morphological redshifts in the case of single-frequency 
interferometers results in a deterioration of the redshift estimation 




real redshift z real 

Figure 14. Redshift determination and error estimation from all three pa- 
rameters a, c and s that followed from wavelet analysis with the sym2- 
wavelet. The upper panel shows the estimated redshift z es t and its error Az 
and the lower panel shows the relative accuracy Az/(1 + z), both as a func- 
tion of redshift z rea i. Here, CMB fluctuations were not taken into account. 
The value of the wavelet moments was set to be q = 3. 
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Table 3. Averaged accuracy of the redshift-determination relative to 1 + z 
based on three parameters derived from the wavelet spectrum of order q 
without the noise contribution from the fluctuating CMB. 



accuracy of a factor close to 1.5, which is caused by imperfections 
of the CMB removal by 5'*-order spline polynomials. 



6 SYSTEMATICS 

SZ clusters would be self-similar and would perfectly follow scal- 
ing relations provided several requirements are fulfilled: ((') virial 
equilibrium (T oc M 2/3 ), (if) structural identity, expressed in equal 
form factors, (Hi) a universal baryon fraction and (iv) the absence 
of heating and cooling processes. Each of these assumptions may 
be challenged and leads to deviations from the self-similar scaling 
relations. While the first two points are included in the numerical 
simulation and are limited by the selection criteria, they increase 
the scatter in the relations between virial quantities, or equally, the 
wavelet parameters a, c and s. Systematic trends of the baryon frac- 
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Figure 15. Redshift determination and error estimation from all three pa- 
rameters a, c and s resulting from wavelet decomposition of the combined 
maps (i.e. with CMB) using the .sym2-wavelet. In the upper panel, the esti- 
mated redshift z es t and its error Az is shown a function of real redshift Zreal- 
In comparison, the relative accuracy Az/(l+z) as a function of z re al is shown 
in the lower panel. Again, the order of the wavelet moments was taken to 
be q = 3. 
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Table 4. Averaged accuracy of the redshift-determination relative to 1 + z 
based on three parameters derived from the wavelet spectrum of order q 
with the noise contribution CMB caused by fluctuations in the CMB. 



tion with cluster mass (see Sect. 6.1) and the formation of cooling 
flows (Sect. 6.2) need to be assessed separately. 

6.1 Influence of tilted scaling relations 

Analyses of X-ray observations carried out by Arnaud & Evrard 
(1999) and Mohr et al. (1999) suggest a weak trend of the clusters 
baryon fraction with cluster mass M and a deviation from the uni- 
versal value f B = Cl b /Q m , which is due to feedback processes like 
galactic winds that more effectively deplete the ICM of baryons 
in low temperature clusters compared to high temperature clusters. 
The variation of the gas fraction, however, is highly model depen- 
dent and reflects the mass and density estimation technique: Es- 
timating masses by means of the virial theorem at fixed density 
contrast, yields a scaling of T cc M 2/3 and requires the baryon frac- 
tion to increase with temperature. Alternatively, using the fit of a 



/3-profile to the cluster, one finds the baryon fraction f B to be nearly 
constant while a steeping of the slope of the M-J-relation is found. 

In either of the above mentioned cases the dependence espe- 
cially of the wavelet parameter a, which is a logarithmic measure 
of the SZ flux T would be increased in more massive clusters and 
would thus increase the scatter in a of a cluster sample at a given 
redshift. The quoted analyses of X-ray data find the baryon fraction 
to show a relative variation amounting to =s 10% at fixed temper- 
ature, i.e. at fixed depth of the potential well for a sample of local 
clusters. Apartfrom the the systematic component, that can in prin- 
ciple beremoved, once high quality X-ray data will improve our 
understanding of this phenomenon and allows proper modelling, 
the stochastic contribution can only be constrained to be at most of 
equal relative influence to AY/T as the scatter in morphology. 

6.2 Cooling flow clusters 

In order to estimate the accuracy of the method outlined above, 
so far we only used adiabatic hydrodynamical simulations which 
lack of cooling processes. Thus we need to address the influence of 
cooling flow clusters on our proposed method. After an analytical 
investigation following Sect. 3.2.1 we compare clusters with and 
without cooling flows and show how the morphological changes in 
cooling flow clusters impacts on the wavelet spectra. 



6.2.1 Analytical wavelet transform of cooling flow clusters 

Instead of a single King profile we assume that the SZ emission of 
a cooling flow cluster can be described by a double King profile for 
reasons of analytical feasibility: 



y(x) =y(r) = Y v, 1 +(- 



(40) 



where the second term describes the additional enhancement owing 
to the cooling flow. Deprojecting this two-dimensional profile by 
means of Pfrommer & EnBlin (2003) yields: 



1 2 2 
p e (R) = - n e (R)kT e (R) = ^ V JL. 

2 cr T 4^ 27IT > 



(l+RVrf) 



3/2 ' 



(41) 



where R denotes the three-dimensional radius and S(a, b) denotes 
the /3-function (Abramowitz & Stegun 1965). Thus we obtain for 
the ratio of the central values of the Comptonisation parameters y t 

^ = ^~i (42) 

yi pi n 2 

where we inserted typical values for cooling flow clusters, p 2 /Pt ~ 
3 and r 2 /ri ~ 1/6. The second order wavelet moment of cooling 
flow clusters can be obtained in analogy to the non-cooling flow 
case: 

Xf{a) = 2nJ dick 5 exp(-fcV) |yir 2 AT (^n) + y 2 r 2 2 K Q (kr 2 f .(43) 

This second order wavelet moment shows an increasing am- 
plitude and a decreasing cutoff parameter compared to the one with- 
out a cooling flow. 



6.2.2 Numerical analysis 

In order to scrutinise these findings we apply our method to adia- 
batically simulated clusters to which we add an enhanced emission 
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Figure 16. The influence of cooling flows to the spectrum of wavelet mo- 
ments, together with the fitting formula (33) for increasing wavelet moment 
order q: q = 2 (squares), q = 3 (circles, solid), and q = 4 (diamonds, 
dashed) for a single cluster without instrumental smoothing. Open symbols 
are values derived from the simulated non-cooling flow cluster, whereas 
filled symbols denote the corresponding cooling flow cluster. 

to mimic the SZ emission of the cooling flow. In Fig. 16 the result- 
ing spectra of wavelet moments are shown together with the fitting 
formula eqn. (33) for increasing wavelet moment order q. 

It can clearly be seen in Fig. 16 that the enhanced emission due 
to the cooling flow yields a slightly higher amplitude of the wavelet 
spectrum on small scales. Extracting information from the wavelet 
spectrum by means of eqn. (33) reveals slightly higher values for 
the amplitude a and smaller values for the cutoff c on the percent 
level. However, this influence is minimised when considering fi- 
nite instrumental resolution particularly for high redshift clusters. 
In any case, if a prominent cooling flow is sufficiently well resolved 
it could be masked and replaced by an interpolation in between the 
mask boundaries. 



7 REDSHIFT ESTIMATION IN A NUTSHELL 

This section shall provide a short summary of how to apply our 
method to an SZ survey for estimating redshifts provided a temper- 
ature map of a patch on the sky with resolved images of clusters. 

• Once a cluster candidate has been localised at a particular po- 
sition of the map this cluster and its ambient field has to be cut out. 
If the number of grid points amounts below 64 2 sampling points, 
the mesh should be refined by interpolation in order to reach dy- 
namical range of approximately two decades. This is important in 
order to provide a sufficiently broad range of scales to be probed by 
the wavelet decomposition. 

• The wavelet spectrum of the map is obtained by wavelet 
transforming the map preferentially using the symlet basis func- 
tions (compare Sect. 3.1). The morphological information con- 
tained within the wavelet spectrum can be extracted by means of 
the model function of eqn. (33). In the case of single-frequency ob- 
servations the ambient CMB field cannot be separated from the SZ 
signal of a cluster. The method described in Sect. 5.8 might be ap- 
plied in order to reconstruct the wavelet spectrum of the pure SZ 
cluster signal. 



• The redshift dependence of the wavelet parameters (amplitude 
a, cutoff c, and slope s) follows the functional form of eqn. (37). 
However, the single model parameters depend on the definitions of 
the particular wavelets and the details of the survey, including dif- 
ferent sources of noise and the cluster detection criteria. The most 
promising way of determining the parameters of the gauge func- 
tions laid down in eqn. (37) would be to derive them from a learn- 
ing set of clusters with known (photometric) redshifts. The final 
redshift estimate of the cluster is most conveniently determined by 
means of maximum likelihood analysis, as described by eqn. (39). 



8 SUMMARY 

In this paper, a method of estimating the redshift of a cluster based 
on the wavelet decomposition of its resolved SZ morphology is pre- 
sented. From a fit to the spectrum of wavelet moments three spec- 
tral parameters are derived, that in turn are non-degenerate and in- 
dicative of cluster distance. These parameters are utilized, through 
a maximum likelihood technique, for estimating the cluster's red- 
shift. In the maximum likelihood technique, empirical gauge func- 
tions describing the wavelet parameter's z-dependence are used. 

First, the method was tested on a simple analytical case: The 
spectrum of Mexican-hat wavelet moments can be derived analyti- 
cally for a King-profile, which is known to describe the Compton-y 
amplitude of clusters well. The spectrum of wavelet coefficients as 
a function of wavelet scale <x, exhibits a break at the cluster scale r c 
and may thus serve as a measure of the cluster's size. Additionally, 
the asymptotic behaviour of the wavelet spectrum in the limit of 
cr » r c and tr <<c r c can be understood. The derivation of wavelet 
moments of order q = 2 is analogous to considering the Fourier 
power spectrum of the Compton-y map, filtered with Fourier trans- 
formed wavelet. The shape of the spectrum of wavelet moments of 
order q = 2 from the analytic calculation is consistent with one ob- 
tained by applying wavelet decomposition to simulated SZ cluster 
maps. 

The method was then applied to set of numerically simulated 
SZ clusters with 1' (FWHM) resolution - comparable to the reso- 
lution of future SZ experiments. The sample comprises 690 cluster 
maps distributed in 23 redshift bins, which is a comparably large 
cluster sample. The clusters are chosen such that they are not in a 
merging state and their SZ image is not too elongated, two crite- 
ria that favour clusters close to virialisation. Additionally, in order 
to simulate single-frequency observations, the cluster maps were 
combined with realisations of the CMB that constitute the main 
source of noise. 

The method was tested for a range of wavelet functions (e.g., 
symlet, coiflet, Daubechies, biorthogonal). The average attainable 
accuracy in estimating redshifts is found to be almost indepen- 
dent of the specific functional form used, although the symlet basis 
yielded the best results. However, the method could benefit from 
improvements concerning the choice of the wavelet basis. For in- 
stance, one could try to construct an optimised wavelet specifically 
for y3-profiles, that yields maximised wavelet coefficients x(n, cr). 

As expected, there is only a weak change in accuracy with 
respect to the order q of the chosen wavelet moment X q (cr). This, 
however, is most likely to change when applying the wavelet anal- 
ysis to noisy images, because for increasing choices of q, uncorre- 
cted noise is suppressed relative to the cluster's signal and concen- 
trating on higher values for q should provide a more robust mea- 
surement of the set of structural parameters a, c and a-. The incre- 
ment of q itself is limited by numerics - this is the case when the 
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wavelet moment X q (cr) is dominated by the largest wavelet expan- 
sion coefficient ^(//, cr), and does not reflect anymore the depen- 
dence on the wavelet scale a. In this limit, the wavelet spectrum 
would exhibit a generic power law behaviour: X q (a) oc a 7 ^ for 
large q. The structural parameters a, c and s were found to depend 
on redshift z by a simple exponential (eqn. (37)). The free param- 
eters in this equation can be determined from a (relatively small) 
sample of SZ cluster images with known redshift. 

The accuracy of determining cluster distances has been as- 
sessed by maximum likelihood estimation. The method yields ac- 
curacies of 4 - 5% relative to 1 + z, which is competitive with pho- 
tometric redshifts, but reaches out to larger distances. At redshifts 
exceeding zi 1, the accuracy is expected to degrade because the 
angular diameter distance ^a(z) starts to level off and thus sets the 
limit of applicability. For single frequency data, the CMB fluctua- 
tions can be removed with a simple polynomial reconstruction ap- 
proach; the accuracy in the redshift estimation is then decreased to 
6 - 7%. 

In this work we have considered two major systematic effects 
that might degrade the accuracy of the method. The first is the vary- 
ing baryon fraction with cluster mass, which has been studied only 
for local cluster samples. While the systematic trend could in prin- 
ciple be corrected for, the stochastic contribution will always add to 
the uncertainty of the distance determination. Another systematic is 
the influence of cooling flows at the cluster's centre. In this case we 
have been able to show that the uncertainty it adds to the redshift 
estimate is very small, mainly because the volume occupied by the 
cooling flow region is limited to the cluster's core. 

Although the result in the distance estimation is stated in terms 
of redshift, it should be emphasised that a specific cosmology is as- 
sumed, which is needed for converting the observables, namely the 
wavelet parameters, to a distance estimate. The distances following 
from the analysis have been expressed as redshifts because of their 
elementary interpretation, but the implicit assumption of an under- 
lying cosmology should be kept in mind when comparing to e.g. 
photometric redshifts. For that reason, the precision of the method 
presented is limited by the accuracy to which the cosmological pa- 
rameters are known. Apart from being a distance indicator, the red- 
shift also plays the role of an evolutionary parameter. 

Comparing this work to the pioneering paper by Diego et al. 
(2003), our expectations concerning the accuracy of morphologi- 
cal redshifts are even more optimistic: Without fitting y6-profiles to 
the observational data, it is possible to describe the cluster's SZ 
morphology by solely relying on wavelet decomposition. Also, we 
describe the spectrum of wavelet moments with a small set of struc- 
tural parameters, that have a lucid physical interpretation, provide a 
non-degenerate distance measurement and enable redshift determi- 
nation owing to their monotonic decline with redshift. Furthermore, 
the redshift dependence of the structural parameters is calibrated 
with the data set itself without relying on prior and simplifying 
assumptions. In spite of the small number of observables consid- 
ered here, the accuracy in the redshift estimation of this method is 
doubled, in comparison with Diego et al. (2003), even for a single 
frequency experiment. 



ACKNOWLEDGEMENTS 

The authors would like to thank Volker Springel and Lars Hern- 
quist who kindly provided us with their numerical simulations and 
Antonaldo Diaferio, Matthias Bartelmann and Simon D.M. White 
for clarifying discussions and many useful comments. 



REFERENCES 

Abramowitz, M. & Stegun, I. A. 1965, Handbook of mathematical 
functions (Dover Books on Advanced Mathematics, New York: 
Dover) 

Arnaud, M. & Evrard, A. E. 1999, MNRAS, 305, 631 
Bartelmann, M. 2001, A&A, 370, 754 

Bennett, C. L., Halpern, M., Hinshaw, G., et al. 2003, ApJ, sub- 
mitted, astro-ph/0302207 
Birkinshaw, M. 1999, Phys. Rep., 310, 98 

Carlstrom, J. E., Holder, G. P., Reese, E. D. 2002, ARA&A, 40, 
643 

Cavaliere, A. & Fusco-Femiano, R. 1978, A&A, 70, 677 

da Silva, A. C, Kay, S. T, Liddle, A. R., Thomas, P. A. 2003, 

astro-ph/0308074 
da Silva, A. C, Kay, S. T, Liddle, A. R., et al. 2001, ApJL, 561, 

L15 

Daubechies, I. & Bates, B. J. 1993, Acoustical Society of America 

Journal, 93, 1671 
Deeming, T. J. 1964, MNRAS, 127, 493 

Diego, J. M., Mohr, J., Silk, J., Bryan, G. 2003, MNRAS, 341, 
599 

Eke, V. R., Cole, S., Frenk, C. S., Patrick Henry, J. 1998, MNRAS, 
298, 1145 

Gradshteyn, I. S. & Ryzhik, I. M. 1994, Table of integrals, se- 
ries and products (New York: Academic Press, — cl994, 5th 
ed. completely reset, edited by Jeffrey, Alan) 
Haiman, Z., Mohr, J. J., Holder, G. P. 2001, ApJ, 553, 545 
Halverson, N. W., Leitch, E. M., Pryke, C, et al. 2002, ApJ, 568, 
38 

Henry, J. P. 2000, ApJ, 534, 565 

Hu, W. 2000, Phys. Rev. D, 62, 43007 

Mohr, J. J., Mathiesen, B., Evrard, A. E. 1999, ApJ, 517, 627 
Monaghan, J. J. & Lattanzio, J. C. 1985, A&A, 149, 135 
Muzy, J. F, Bacry, E., Arneodo, A. 1993, PRE, 47, 875 
Navarro, J. F, Frenk, C. S., White, S. D. M. 1995, MNRAS, 275, 
720 

Pfrommer, C. & EnBlin, T. 2003, A&A, in press, astro- 

ph/0306257 
Rephaeli, Y. 1995, ARA&A, 33, 541 
Seljak, U. & Zaldarriaga, M. 1996, ApJ, 469, 437 
Silk, J. 1968, ApJ, 151,459 

Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJ, submitted, 

astro-ph/0302209 
Springel, V. & Hernquist, L. 2002, MNRAS, 333, 649 
Sunyaev, R. A. & Zel'dovich, I. B. 1972, Comments Astrophys. 

Space Phys., 4, 173 
— . 1980, ARA&A, 18, 537 
van den Bosch, F. C. 2002, MNRAS, 331, 98 
Wechsler, R. H., Bullock, J. S., Primack, J. R., Kravtsov, A. V., 

Dekel, A. 2002, ApJ, 568, 52 
White, M., Hernquist, L., Springel, V. 2002, ApJ, 579, 16 
Zaroubi, S., Squires, G, de Gasperis, G, et al. 2001, ApJ, 561, 

600 

Zaroubi, S., Squires, G, Hoffman, Y, Silk, J. 1998, ApJL, 500, 
L87+ 

Zhao, D. H., Mo, H. J., Jing, Y. P., Borner, G. 2003, MNRAS, 339, 
12 

This paper has been typeset from a TgX/ ETj-X file prepared by the 
author. 



